Signatures of Electronic Nematic Phase at Isotropic-Nematic Phase Transition 
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The electronic nematic phase occurs when the point-group symmetry of the lattice structure is 
broken, due to electron-electron interactions. We study a model for the nematic phase on a square 
lattice with emphasis on the phase transition between isotropic and nematic phases within mean 
field theory. We find the transition to be first order, with dramatic changes in the Fermi surface 
topology accompanying the transition. Furthermore, we study the conductivity tensor and Hall 
constant as probes of the nematic phase and its transition. The relevance of our findings to Hall 
resistivity experiments in the high-T c cuprates is discussed. 

PACS numbers: 71.10.Hf,72.15.Eb,72.80.Ga 



X 



I. INTRODUCTION 

Recently, novel spatially inhomogeneous and/or 
anisotropic phases have been found to occur in various 
condensed matter systems. Similar to classical liquid 
crystals, lj these quantum phases can be classified ac- 
cording to their broken symmetries. One such phase is 
the smectic phase, often referred to as a stripe phase. 
This phase can be viewed as an unidirectional charge den- 
sity wave phase which spontaneously breaks translational 
symmetry along one direction. Experimental evidence 
for smectic phases has been found in several transition- 
metal oxide compounds. In particular, a stripe crystal 
phase has been observed in manganese oxide compounds 
by electron diffraction. |2( Moreover, neutron scattering 
experiments have provided indirect evidence for stripe 
phases in the high-T c cuprate superconductors. 

Another such phase is the nematic phase - - an 
anisotropic phase with broken rotational symmetry, but 
with translational symmetry in tact. This phase can 
be thought of as a smectic phase "melted" via the pro- 
liferation of topological defects. Nematic phases have 
been argued to be the source of the anisotropic trans- 
port observed in quantum Hall systems in high Landau 
levels 0,01, and to be relevant to the high-T c cuprate su- 
perconductors as well.[6j However, an understanding of 
the transition from the smectic to nematic phase due to 
the proliferation of the topological defects still remains 
to be achieved. 

Recently, a complementary approach to studying the 
nematic phase was suggested. 0,1^ In this approach, the 
nematic phase is obtained from an isotropic liquid by 
breaking rotational symmetry. In the continuum limit, 
it was shown that the Goldstone mode resulting from 
breaking rotational symmetry leads to non-Fermi liquid 
behavior in the single-particle lifetime. Exotic super- 
conducting instabilities mediated by the collective mode 
in the nematic phase were also studied. [|| 

While the consequences of the nematic order in the 
recent complementary approach were studied in Ref. 
ITol ITTl we still lack an understanding of the phase tran- 
sition between the isotropic and nematic phases. In this 



work, we study a model for the nematic phase on a square 
lattice, with emphasis on the phase transitions between 
the isotopic and nematic phases. Furthermore, we study 
the conductivity tensor and the Hall constant as probes 
of nematic phase and its transition. 



II. MODEL HAMILTONIAN AND PHASE 
TRANSITIONS 



The Hamiltonian we consider is 
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In Eq. ^ £k is the single-particle dispersion; F 2 (q) is the 
inter-electron interaction strength, which can be written 
as 
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Q(q) is the quadrupole density tensor 
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where N is the number of lattice sites. In what follows, 
we take £k to arise from a tight-binding model on a square 
lattice with nearest-neighbor (i) and next-nearest neigh- 
bor (t') hopping integrals 



£k = — 2t(cos k x + cos k y ) — At' cos k x cos k y 



where fj, is the chemical potential, and we have set the 
lattice spacing to unity, a = 1 . 

Within a mean field approximation, we have the self- 
consistency equations 
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FIG. 1: The nematic order A as a function of chemical poten- 
tial, fi for \F 2 \/(2t) = 0.55 and if = -0.4i (2t = 1) at T = 0. 
/i = —0.4 corresponds to the half filling, x — 

where /(£k) is the Fermi function, and n is the electron 
density. The expectation values (Q xx ) and (Q X y) lead to 
a modification of the single-particle dispersion. Defining 

X = \F 2 \(Q xx )/{2t) and A' = \F 2 \(Q xy )/(4t') , (3) 

the single particle dispersion becomes 

£k = — 2t [(1 + A) cos k x + (1 — A) cos k y ] 

— 4i'(cos fca; cos ky + A' sin k x sin — . (4) 

For given values of fi, F 2) and i', we solved for A, A' and 
the density n = 1 — x (x is the hole concentration) using 
Eq. [3 For all values of the parameters we considered, we 
found A' = 0. 

In Fig. ^ we show the behavior of the nematic order 
A as a function of chemical potential, /i for F 2 /(2t) = 
-0.55, if = -OAt, and 2t = 1 at T = 0. An important 
feature of Fig. is the sharp onset and offset of nematic 
order at critical values of fi: fi c — —0.57 and —0.94. The 

free energy near the transition 0.61 < fi < —0.53 

and —0.96 < fi < —0.89 — shows three minima at finite 
±A and A = 0. This indicates a region of coexistence 
of the isotropic and nematic phases: it reflects the nega- 
tive quartic term in the Ginzburg-Landau free energy. |12j 
Hence, the transition between the isotropic and nematic 
phases is first order. The electron density as a function 
of [i shows a discontinuity at these critical values of ^t c , 
which is a signature of the first order transition. 

To better understand the behavior of the nematic or- 
der, as well as the robustness of our result, we considered 
the dependence of A on the parameters F 2 and t'. Upon 
increasing \F 2 \, we found the maximum value of A to 
increase. Also, increasing \F 2 \ shifts the critical values 
of /i, making the region over which nematic order exists 
wider. Furthermore, we found an absence of nematic or- 
der for all doping concentrations if \F 2 \/(2t) < 6 x 10~ 3 
(within the numerical precision) for if = —OAt. By vary- 
ing t' , we found the structure shown in Fig. ^to shift. 



In particular, when t' = 0, A is maximal at /i = (i.e. 
half-filling) and is a symmetric function of /x. This is a 
consequence of the particle-hole symmetry present when 
t' = 0. Though the critical values of /i and the maximal 
value of A depend on the value of F 2 and t' , we would 
like to stress that the qualitative features of the phase 
transition between the isotropic and nematic liquids are 
not sensitive to the choice of parameters. In what fol- 
lows, we will present results for F 2 /(2t) = —0.55, and 
if = — OAt (2t = 1). For small hole concentrations, this 
gives a hole-like Fermi surface centered about (tt, tt) (see 
Inset (a) of Fig. EJ). Such a Fermi surface is consistent 
with angle resolved photoemission spectroscopy^ anc l 
Hall constant measurements 0] in the high-T e cuprates. 

III. CONDUCTIVITY TENSOR AND THE 
HALL CONSTANT 

Above, we saw that the nematic phase occurs in a fi- 
nite regime dopings. Furthermore, we found a dramatic 
change in the Fermi surface topology upon developing 
nematic order. (See the Inset of Fig- El) I n this section, 
we study the conductivity tensor to probe the nematic 
order, as the conductivity tensor, especially the hall con- 
ductivity is sensitive to the topology of the Fermi surface. 

Within a Boltzmann equation approach, the longitu- 
dinal conductivities are (in units where % = 1) 

where i = x or y, e is the electron's charge, Vi is the 
component of the Fermi velocity along the i-axis, and Tk 
is the transport lifetime. In the presence of a magnetic 
field H parallel to the £-axis, the Hall conductivity is[l5L 

m 

a xy = 2e 3 ff— f J K r k) ( v k x V k ) 2 (v x T k ) ■ 

In what follows, we will report results for the case Tk = 
t = const. Some discussion regarding the form of Tk is 
given in the following Section. 

Using the single-particle dispersion in Eq. 01 with A 
determined from Eq. |2 we plot the conductivity tensor 
vs. [i in Fig. For comparison, the conductivity tensor 
without nematic order (A = 0) is shown as a dashed line. 
When nematic order develops, the Fermi surface topol- 
ogy changes dramatically. More specifically, the hole- 
like Fermi surface with nematic order becomes an open 
Fermi surface upon developing nematic order. In the ne- 
matic state, the Fermi surface is rather flat along the 
y-direction. (See Inset (b) of Fig. |3) Hence, the hole's 
motion along the ^-direction is enhanced, while it is sup- 
pressed along the y-direction. These changes in the Fermi 
surface are reflected in the longitudinal conductivities, 
Figs. Ufa) and |2{b) — a xx increases suddenly with the 
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FIG. 2: The conductivity tensor as a function of chemical 
potential fi for \F 2 \/{2t) = 0.55 and t' = -OAt (2t = 1). 
(a) a xx vs. jj,, (b) <jyy vs. /i, and (c) the Hall conductivity 
a xy vs. fx. a xx and a^j, are in units of 2e 2 (2i) 2 ; a xy is in 
units of 2e 3 H(2t) 3 . The dramatic change is due to the phase 
transition to nematic state. The dashed lines indicate the 
behavior of conductivities without the nematic order. 



onset of nematic order (at fi = —0.57), while a yy shows a 
sudden decrease. The Hall conductivity, Fig.^c), shows 
a dramatic drop at /i = —0.57 as well. However, its 
change is small compared with a xx and a yy . This occurs 
because the Hall conductivity is not only sensitive to the 
Fermi surface's topology, but also its curvature. As 
shown by the dashed line in Fig. [3f c) , the Hall conduc- 
tivity already decreases fairly rapidly without nematic 
order. The dotted lines at the transition indicate that 
Boltzmann equation, Eq. [S] is not applicable near the 
transition, since there are two phases with different den- 
sities coexist. 

Another quantity of interest is the Hall constant nn, 
which measures the number of charge carriers, and is 
related to the Hall resistivity. For weak magnetic fields, 
the Hall resistivity is 
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which defines tir in terms of p xy . For the case of a single- 



FIG. 3: Hall constant as a function of chemical potential, fi 
for F 2 /(2t) = -0.55 and t' = -OAt (2t = 1). The dashed line 
is the Hall constant without the nematic order. The inset 
shows the the Fermi surface for the points labeled by arrows 
in the main figure. 



band Fermi liquid model with a circular (or elliptical) 
Fermi surface, nn is the density of charge carriers per 
unit cell i.e. the area of the Fermi surface. 

In Fig. we plot the Hall constant, n# as a function 
of chemical potential, /i. For comparison, the Hall num- 
ber without nematic order (A = 0) is shown as a dashed 
line. Notice the sharp drop of nn at fi = —0.57, where 
nematic order becomes finite. It is expected that nn is 
almost constant in the region where nematic order ex- 
ists, —0.94 < fi < —0.57. This can be understood from 
the conductivity tensor in Fig. [21 — the decrease of hole 
motion along the y direction is compensated by the in- 
crease of hole motion along x direction, giving a product 
&xx&yy I &xy ~ ^ff which is nearly constant. 



IV. DISCUSSION 

To begin with, our results for the phase transition 
in Fig. n were obtained within a mean-field approxima- 
tion. Once fluctuations about the mean-field state are 
included, we expect the range over which the nematic 
phase exists to be reduced. Furthermore, fluctuations 
will reduce the value of A. The effect of fluctuations on 
the nature of phase transition will be addressed in near 
future. Another approximation in our treatment was the 
neglect of collective modes. In the continuum limit of 
Eq. ^ the breaking of rotational symmetry leads to a 
Goldstone mode; scattering from this Goldstone mode 
gives rise to non- Fermi liquid behavior. On a lattice, 
however, this collective mode is gapped; we expect scat- 
tering from the collective mode to have little effect on 
the dc conductivities. 

The subject of disorder deserves some discussion. 
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Though disorder was introduced to understand the be- 
havior of the conductivity through the transitions, we 
did not address its effects on the transition between the 
isotropic and nematic phases, nor did we address its con- 
sequences on the nematic phase itself. 0| I n the absence 
of disorder, the mean-field treatment of the Hamiltonian 
Eq. n gave the transition between the isotropic and ne- 
matic phases to be first order. For systems that undergo a 
first order transition without disorder, it was shown that 
the system with quenched disorder could be mapped onto 
a random field Ising model. [18| Based on this mapping, 
it was shown that in d < 2 dimensions, the coexistence 
at the phase transition disappears for arbitrarily small 
amount of disorder; in d > 2 dimensions this occurs 
by the application of a finite amount of disorder. For 
the transition in this work, however, there are gapless 
fermions presents on both sides of the transition. There- 
fore it is unclear to us how much the results obtained 
from the mapping may apply. We hope to address these 
issues in the future. 

It is also important to note that the conductivity tensor 
depends of the form of the transport scattering rate, and 
hence on the disorder. Realistic models for disorder usu- 
ally give transport scattering rates which are momentum 
dependent. In particular, in most high-T c materials, dis- 
order is mainly due to the dopant ions, which lie outside 
of the CuOi planes. As a result, these dopant ions cause 
primarily forward scattering. In the above discussion, we 
described our results for the case Tk = r. However, we 
have also considered the case where Tk ~ |vk|. Such 
a form arises for out-of-plane disorder, causing primarily 
forward scattering. 0] When A = 0, we found the overall 
features of nu to be very different from Tk = r. However, 
the effect of the nematic order on the Hall constant gives 
the same result as when Tk = r — there is a sharp drop 
in the Hall constant in the nematic phase compared with 
the isotropic phase. 

Our findings for the Hall constant have poten- 
tial relevance to the experiments of the normal 
state Hall resistivity in the high-T c superconductor, 
Bi 2 Sr 2 _2 : La a; Cu06+5.[23, Hl| In these experiments, the 
Hall constant was measured in the normal state at low 
temperatures, by suppressing superconductivity with a 



large pulsed magnetic field of 60 T. They observed satu- 
ration of the Hall constant at low temperatures, suggest- 
ing conventional transport properties. However, the Hall 
constant as a function of doping concentration shows a 
maximum at x — 0.15, and decreased sharply in under- 
doped regime. In overdoped regions, there is a regime 
near optimal doping x ~ 0.16, where the Hall constant 
shows the dip — the decrease is about 20% (~ 0.2) of its 
maximum value (~ 1.0) at x = 0.15. The dip of the Hall 
constant around x = 0.16 becomes more pronounced at 
low temperatures (< IK). This suggests that the ne- 
matic phase may be a candidate to explain the anomaly 
observed in the Hall constant in high T c cuprates. 

In conclusion, to study the nematic phase we consid- 
ered the Hamiltonian in Eq.Uwith quadrupolar density- 
density interactions on a square lattice. Within a mean- 
field analysis, we found the nematic phase to occur in 
a finite region of dopings. Moreover, the transition be- 
tween the isotropic and nematic phases was found to be 
first order; a dramatic change in the Fermi surface topol- 
ogy accompanied this transition. As the conductivity is 
sensitive to the topology of the Fermi surface, we investi- 
gated the effects of nematic order on transport. We found 
the conductivity tensor and Hall constant to show sig- 
nificant changes at the transitions between the isotropic 
and nematic phases. If observed, such signatures would 
give evidence for nematic order; it would give evidence 
for the approach pursued in this paper as an appropri- 
ate description of nematic phases in physical systems. 
Finally, we discussed the effect of disorder and the scat- 
tering from collective modes on the phase transitions and 
on the transport properties at the transitions. 

Acknowledgments We would like to thank S. A. Kivel- 
son, V. Oganesyan P. Coleman, Y-B Kim, and E. 
Fradkin for insightful discussions, and Yoichi Ando 
for useful comments. We also thank P. Fournier 
and W. Metzner for pointing out Ref. [lj and Ref. 
0, respectively. This work was supported by the 
NSERC of Canada(HYK,EHK,CHC), Canada Research 
Chair(HYK), Canadian Institute of Advanced Re- 
search(HYK), Alfred P. Sloan Research Fellowship 
(HYK), and Emerging Material Knowledge program 
funded by Materials and Manufacturing Ontario (HYK). 



[1] Physics of Liquid Crystals, P. de Gennes ( Oxford Uni- 
versity Press, 1987). 

[2] S. Mori, C. H. Chen, and S. W. Cheong, Nature, 392, 
473 (1998). 

[3] J. Tranquada, B. J. Sternlieb, J. D. Axe, Y. Nakamura, 

and S. Uchida, Nature, 375, 561 (1995). 
[4] M. P. Lilly, K. B. Cooper, J. P. Eisenstein, L. N. Pfeiffer, 

and K. W. West, Phys. Rev. Lett. 82, 394 (1999); R. 

R. Du, D. C. Tsui, H. L. Stormer, L. N. Pfeiffer, K. 

W. Baldwin, K. W. West, Solid State Comm. 109, 389 

(1999); K. B. Cooper, MP. Lilly, J. P. Eisenstein, L. N. 

Pfeiffer, K. W. West, Phys. Rev. B 65, 241313 (2002). 



[5] E. Fradkin and S. A. Kivelson, Phys. Rev. B 59, 8065 
(1999). 

[6] S. A. Kivelson, E. Fradkin, and V. J. Emery, Nature 393, 
550 (1998). 

[7] V. Oganesyan. S. Kivelson, and E. Fradkin, Phys. Rev. 

B 64, 195109 (2001). 
[8] C. J. Halboth, and W. Metzner, Phys. Rev. Lett. 85, 

5162 (2000); A. Neumayr, and W. Metzner, Phys. Rev. 

B 67, 035112 (2003); V. Hankevych, and F. Wegner, Eur. 

Phys. J. B 31, 333 (2003). 

[9] Y. B. Kim and H.-Y. Kee, cond-mat/0204037 
[10] H.-Y. Kee, Phys. Rev. B, 67, 073105 (2003). 



5 



[11] S. A. Kivelson, E. Fradkin, V. Oganesyan, I. P. Bind- 
loss, J. M. Tranquada, A. Kapitulnik, and C. Howald, 
cond-mat/0210683 to be published in Rev. Mod. Phys. 

[12] C. H. Chung, H.-Y. Kee, and E. H. Kim, unpublished. 

[13] See the review article: J. C. Campuzano, M. R. Nor- 
man, and M. Randeria, in in Physics of Conventional and 
Unconventional Superconductors, edited by K. H. Benne- 
mann and J. B. Ketterson ( Springer- Verlag). 

[14] H. Y. Hwang, B. Batlogg, H. Takagi, H. L. Kao, J. Kwo, 
R. J. Cava, J. J. Krajewski, and W. F. Peck, Jr., Phys. 
Rev. Lett. 72, 2636 (1994). 

[15] N. P. Ong, Phys. Rev. B 43, 193 (1991). 

[16] J. Ziman, Electrons and Phonons (Clarendon, Oxford 
1960). 

[17] The effect of the disorder on electronic liquid crystal is 



one of the interesting issues. For the classical liquid crys- 
tal in the continuum, see L. Radzihovsky and J. Toner, 
Phys. Rev. B 60, 206 (1999) 

[18] J. Cardy, cond-mat/9806355 See the references therein. 

[19] H.-Y. Kee, Phys. Rev. B 64, 012506 (2001). 

[20] F. F. Balakire, J. B. Betts, S. Ono, T. Murayama, Y. 
Ando, and G. S. Boebinger, in Physical Phenomena at 
High Magnetic Fields-IV edited by G. Boebinger, Z. Fisk, 
L. P. Gorkov, A. Lacerda, and J. R. Schrieffer (world 
Scientific, Singapore, 2002). 

[21] A theoretical interpretation was given in S. Chakravarty, 
C. Nayak, S. Tewari, and X. Yang, Phys. Rev. Lett. 89, 
277003 (2002). 



